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Abstract 



For some time, the theoretical result for the transition temperature of a di- 
lute three-dimensional Bose gas in an arbitrarily wide harmonic trap has been 
known to first order in the interaction strength. We extend that result to sec- 
ond order. The first-order result for a gas trapped in a harmonic potential 
can be computed in mean field theory (in contrast to the first order result for 
a uniform gas, which cannot). We show that, at second order, perturbation 
theory suffices for relating the transition temperature to the chemical poten- 
tial at the transition, but the chemical potential is non-perturbative at the 
desired order. The necessary information about the chemical potential can 
be extracted, however, from recent lattice simulations of uniform Bose gases. 

I. INTRODUCTION 

Consider a dilute three-dimensional gas of bosons, all identical, in an external harmonic 
trapping potential 

V(x) = l -m(u 2 x x 2 + u 2 y 2 + u 2 z z 2 ), (1.1) 

where m is the mass of each boson. For this system to have a sharp, well-defined phase 
transition, we need to formally take the infinite volume limit of u x , u y , u z — > while keeping 
the central density of Bose particles finite and non-zero at the transition. As we'll briefly 
review below, the central density n at the transition scales as n ~ A 1//2 /a^ Q , where N is the 
total number of Bose particles in the trap, and where 



h \ 1/2 ( h \ 1/2 / h \ 1/2 
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is the volume scale of the ground-state wave function. (See also ref. |J for a review.) So the 
appropriate infinite volume limit is u x , u y , u z — > with Nuj x uj y uj z held fixed. 

At low energies, the relevant measure of the strength of interactions is the 2-body scat- 
tering length a. We will assume that interactions are repulsive (a > 0). We will study the 
transition temperature T c for Bose-Einstein condensation of a dilute single-species gas as 
a function of the total number of particles N in the trap, in the infinite volume limit just 
discussed. One might naively anticipate there to be an expansion of the form 



T C {N) = T (N) 



a /a x 2 
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where T is the ideal gas result and I is some characteristic length of the ideal gas system. As 
we'll review below, the appropriate length scale for a trapped Bose gas is the typical inter- 
particle separation I ~ n -1 / 3 ~ N^^a^ at the center of the trap. The coefficient c\ of the 
expansion for T c has been known for several years 0. In this paper, we calculate the next 
correction. As we'll discuss, this is the furthest one can go in the expansion without more 
information about interactions than just the scattering length. We'll find that c 2 depends 
logarithmically on a /I: the actual expansion is of the form 



T C (N) = T (N) 



a- 2 
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and we shall calculate the constants c' 2 and c 2 '. (The appearance of a related logarithm for 
uniform gases has been qualitatively discussed in Ref. ||. For a calculation of the second- 
order relationship between T c and the central density n in an arbitrarily wide trap, which is 
also the relationship T c (n) for a uniform gas, see Ref. Q.) 

Some aspects of the Bose-Einstein condensation phase transitions are perturbatively cal- 
culable, and others are not. In a dilute Bose gas, the physics of fluctuations associated with 
relatively short distance scales is perturbative, while that associated with critical behavior 
on relatively long distance scales is not. In the case of a uniform Bose gas (that is, a Bose 
gas in an infinite square well potential rather than a harmonic potential), the first-order shift 
in T c is sensitive to critical fluctuations and so is non-perturbative. That shift has recently 
been calculated using lattice simulations and has previously been estimated in a wide 
variety of ways H8Hl4|. In contrast, the first-order shift for a gas trapped in a harmonic 
potential (parametrized by Ci) is calculable using perturbation theory 0. As we shall see, 
the second-order logarithmic coefficient c' 2 is also calculable in perturbation theory, but the 
constant c 2 under the log is not. We shall calculate c 2 by relating it to measurements that 
have been made in lattice simulations of the phase transition in three-dimensional 0(2) field 
theory j7|. 

We should emphasize that expansions of physical quantities in a/ 1 cease to correspond 
to perturbative expansions in a/ 1, once one reaches the orders we have asserted are non- 
perturbative. The failure of perturbation theory in describing generic second-order tran- 
sitions has been known for decades. This breakdown typically manifests in perturbation 
theory as the appearance of infrared infinities in the coefficients of the perturbative expan- 
sion beyond a certain order. 

There is a simple way to relate the problem of a Bose gas in an arbitrarily wide harmonic 
potential with that of a uniform Bose gas. In the infinite volume limit ui x , to y , ui z — > of 
the harmonic trap problem, the trapping potential becomes everywhere locally flat over any 
fixed distance scale (such as the typical inter-particle spacing). Locally, the problem can 
then be treated as a uniform gas in the presence of a x-independent potential, and an x- 
independent potential can be absorbed into a redefinition of the chemical potential. For 
example, if the original chemical potential was /2, then the effective chemical potential at 
a position x is fx — V(x). For arbitrarily wide traps, the total number of particles in the 
system is then related to chemical potential and temperature by 
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N = Jd 3 xn(T, Jl- V(x)), (1.5) 

where n(T, /i) is the uniform gas result for the number density at a chemical potential fi. 

In a trap, the effective chemical potential p, — V(x) is highest at the center, where 
V(x) = 0, and this is where the condensate first forms as the system is cooledj] For a 
uniform gas, let /2 C (T) be the critical value \x of the chemical potential at a given temperature 
T . Then (|1.5|) becomes 

N = J An(T c , fi c (T c ) - V(x)). (1.6) 

If we knew n(T, //) and fi c (T) for a uniform gas, we could then use (|1.6|) to solve for T c for 
a gas of N particles in an arbitrarily wide trap. 

In the next section, we review in more detail the physical scales of the problem and 
explain why, for the purposes of applying (|1.6|) to second order, it is adequate to use pertur- 
bation theory for the uniform gas result n(T,[i). We also explain why perturbation theory 
is inadequate to find the uniform gas result /2 C (T) at second order. The second-order pertur- 
bative result for n(T, //) can be extracted from the literature JTSyiB"! , and in Section pTT| we 



step through the simple exercise of applying that old result to determine the relation (|1.6| 



between N, p, c , and T c at second order. Then, in Section [TV], we take on the less trivial step 
of showing how the second-order value of pL c (T) can be related to existing results from lattice 
simulations of 0(2) scalar field theory in three dimensions. We put everything together in 
Section |V|, giving our final answer for the second-order term of the expansion ( |1.4| ) of T c . 
In Section |VI], we discuss the nature of yet higher-order corrections and explain why they 
require more knowledge of 2-body scattering than just the scattering length a. In Section 
[VI l[ we briefly discuss parametrically how wide a trap must be for our "arbitrarily wide 



trap limit" results to be valid at second order. Finally, we conclude in Section |VIII with 



a brief example of how big the second-order effects might be in a particular experimental 
situation. Various details and diversions are saved for appendices, including a modern field- 
theory rederivation and verification of the old perturbative result for n(T, /i) that we take 
from Huang, Yang, and Luttinger 0,|16| . 



II. SCALES AND EFFECTIVE THEORIES 
A. The uniform gas 

Before proceeding to a Bose gas in a harmonic trapping potential, let's first review the 
basic scales relevant to the phase transition of a uniform gas. The generic condition that 
the gas is dilute is that the two-particle scattering length a be small compared to the typical 
inter-particle separation I ~ n^ 1 ^ 3 , where n is the number density. The Bose-Einstein 
condensation phase transition occurs when the typical de Broglie wavelength 



1 A reminder about signs: Recall that, for a uniform Bose gas, /i is negative at high temperature and 
increases (moves towards zero) as the system is cooled. 
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A = h^2ir(3/m, (2.1) 

becomes of order the inter-particle separation I. Then a < A ~ n^ 1 ^ 3 . 

At the phase transition, the interaction can be treated perturbatively for analyzing short- 
distance physics but, as with most second-order phase transitions, interactions cannot be 
treated perturbatively for analyzing long-distance physics. A distance scale that will be of 
interest is the dividing line between these two regimes. As we shall review below, this scale 
is A 2 /a ~ n _2//3 /a. At the transition, there is then a hierarchy \ 2 /a ^> A ^> a of physically 
relevant distance scales for a dilute Bose gas. 

We will now briefly review the description of the dilute Bose gas system in terms of 
effective field theories, and we'll then turn to the effective field theory description relevant 



to the long distance physics at the critical point |L2| . This will provide a clean way to review 
the origin of the non-perturbative scale A 2 /a, and we will need to make use of such effective 
theories later in our discussion of the critical chemical potential fi c (T) for a uniform gas. 



B. The action 

It is well known that, at distance scales large compared to the scattering length a, 
an appropriate effective theory for a dilute Bose gas is the second-quantized Schrodinger 
equation, together with a chemical potential fi that couples to particle number density 
and a \ip\ 4 contact interaction that reproduces low-energy scattering |IJ. The corresponding 
Lagrangian is 

C = r (ihd t + f- V 2 + fi - V(x)) i> - (^) 2 - (2.2) 

y 2m J m 

The identification of the coefficient of the (i/j*^) 2 interaction with 2irh 2 a/m is technically only 
valid at leading order in the interaction strength but, as we'll review later, doesn't change at 
second order if one uses dimensional regularization ||17|| . We'll also later discuss (in section 



[VTD the size of corrections to the effective theory due, for instance, to energy dependence 
of the cross-section or 3-body interactions. It will turn out that such corrections can be 
ignored for the purpose of computing T c to second order. 

To study (|2.2| ) at finite temperature, apply the standard imaginary time formalism, 
so that t becomes —it and imaginary time r is periodic with period h/3 = h/k^T. The 
imaginary-time action is then 



s, = F dr I if'.r 



^(hd T -^V 2 -^ + V(x)) il> + (^) 2 
V 2m / m 



(2.3) 



As usual, the field ip can be decomposed into imaginary-time frequency modes with Mat- 
subara frequencies u n = 2im/Tif3. 



C. Non-perturbative physics in the uniform gas 

We'll now specialize the preceding to the uniform g V(x) = and will discuss 

the system at or close to the critical point. For distances large compared to the thermal 
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wavelength ( |2.1|) and sufficiently near the transition so that <C T, the non-zero Matsubara 
frequencies decouple from the dynamics, leaving behind an effective theory of only the zero- 
frequency modes ipo, with the action becoming 



h- 1 S l Jd 3 x 



( h 2 2 \ 2irh 2 a . 
Wq ~7T~ V -PeB)V>0 + (Wo 

\ 2m I m 



(2.4) 



up to corrections that again, as we will discuss later (in section |VlD , do not affect T c at 



second order. Eq. ( |2.4j ) can be interpreted, if desired, as the f3H of a classical 3-dimensional 
field theory. 

Finally, it is convenient to rewrite ipo = <f)\4ft/ A so that the effective action becomes a 
conventionally normalized U(l) field theory of a complex field 4>: 



S = I il'.V 

where 



u 

+ r<p*<p + - 
o 



(2.5) 



96n 2 a 

« = — ■ (2.6) 

We will henceforth refer to this effective theory as the "three-dimensional" effective theory, 
while referring to the original imaginary time theory ( |2.3| ) as the "3+1 dimensional theory" 
(for three space plus one time dimension). By writing <fi = (0i + z</>2)/v2, the three- 
dimensional effective theory may equivalently be interpreted as an 0(2) theory of two real 
fields with interaction (w/4!)(0f + 0|) 2 - The relationship of r to the chemical potential \x and 
other parameters of the original theory is a little more subtle, because the 4>*4> interaction 
is associated with an ultraviolet (UV) divergence of the three-dimensional theory that has 
to be renormalized. We will discuss this relationship in detail when we analyze /2 C (T) in 



section [TV]. For the moment, these details are unimportant. 

There will be a line in the (//, T) plane that corresponds to the Bose-Einstein condensation 
phase transition. In the long-distance effective theory ( |2.5| ), that line will correspond to a 
line in the (r, u) plane. If we think of this line as determining r in terms of u, then the 
only physical scale in the problem of studying this effective theory at the transition is u. By 
dimensional analysis, the distance scale of non-perturbative physics is therefore 1/u ~ A 2 /a, 
as asserted earlier. 

It will be useful to understand how far away from the transition one needs to go, as 
measured by p, c — /i at T = T c , in order for the physics on all scales to be perturbative. This 
will happen when the correlation length £ is small compared to the scale 1/u ~ X 2 /a of non- 
perturbative physics. We can determine this condition on £ with a perturbative analysis. 
In fact, it is sufficient to consider a simple Gaussian (i.e. tree-level) approximation, where 
/i e fr in the effective three-dimensional theory (|2.4|) is naively taken to be fi, corresponding to 
r = —2mfi/h 2 in the rescaled effective theory fl2.5p . In Gaussian approximation, /2 C = 0. The 
correlation length, in Gaussian approximation, is £ ~ r~ 1,/2 , and so the condition ^ C 1/k 
becomes 

Ti 2 u 2 h 2 a 2 

fic ~ fl > 7T ■ 2.7 
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a 


scattering length 


a 


a 


I ~ fi-l/3 ^ A 


inter-particle separation at trap center; 
ineimai waveiengtn 


iV- 1/6 a ho 


I 




smallest wavelength of non-perturbative 
fluctuations near center of trap 


AT— 1/3 2 / 


12 1 

I /a 


a ho 


size of the ground state (condensate) 




nvh 




size of non-perturbative region 




N^ 3 a 


-^cloud 


size of entire gas cloud 


iV 1/6 a ho 





TABLE I. Distance scales for a dilute Bose gas in an arbitrarily wide harmonic trap at the phase 
transition. The scales are ranked in ascending order. Entries should be interpreted as representing orders of 
magnitude (parametric dependence) and not as precise definitions and equalities. The first column gives our 
notation for each scale. The third column shows how the scales depend on the "experimental" parameters 
a, to, N, and lu, where aho = (h/muj) 1 ^ 2 . The last column shows a simple rewriting that makes the ordering 
of scales clear, given that our assumed limits can be phrased as a fixed; a <C / (diluteness); and I fixed with 
N — > oo (arbitrarily wide trap). 



(An equivalent condition was discussed in the original work || on the first-order result for 
T c in a trap.) We'll see later, in our more thorough discussion of the relationship between 
r and fi in section |V[ that corrections to the Gaussian approximation do not change this 
conclusion. 

Finally, note that, by dimensional analysis, the non-perturbative contribution to the 
critical value of r in the three-dimensional 0(2) effective theory ( |2.5| ) must be of order 
u 2 . The Gaussian approximation's identification of r with -2m/i/^ 2 then suggests that the 
non-perturbative contribution to the critical value fl c is of order h 2 u 2 /m, which is second 
order in a. As we'll see in section [Iy|, this conclusion is correct. This is the reason that, in 
order to calculate T c to second order, we must account for non-perturbative physics in the 
determination of Jx c . 



D. Gas in a harmonic trap 

Now we turn to reviewing scales in a harmonic trap. One of the main points of this 
exercise will be to determine the size of the region (at the transition) where the physics is 
non-perturbative, relative to the size of the trapped gas cloud as a whole. This will allow 
us to determine to what order one can use perturbative calculations to relate N, T, and fi 
via (|i~5|). 

For simplicity, we'll assume in this discussion that uj x ~ uj y ~ uj z . The relevant distance 
scales for a dilute Bose gas in an arbitrarily wide harmonic trap, at the transition, are 
summarized in Table | in ascending order. Most of this is just review of simple, standard 
results Jl], except for the scales of non-perturbative physics in a harmonic trap, which we 
haven't seen clearly discussed before. 

First, let's review the size and density of the cloud of Bose particles at the phase transi- 
tion. As we'll reproduce below, most of the particles in the trap are in the classical regime, 
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and we can use the classical equipartition theorem to find the width -R c ioud of the cloud: 
|mco> 2 x 2 ~ hkftT yields -R c ioud ~ (Pmuj 2 )^ 1 / 2 . The central density of particles is then of or- 
der n ~ iV/i? 3 loud ~ N(f3muj 2 ) 3 / 2 , and the separation of particles at the center of the trap is of 
order / ~ n -1 / 3 ~ iV _1//3 (Pmuj 2 )^ 1 ^ 2 . The phase transition occurs when this separation is of 
order the thermal wavelength ( |2.1| ) , giving k B T ~ N^htu, and so / ~ n 1 ^ ~ A ~ N 1//6 Oh , 



as claimed in Table |. The fact that k B T ~ N^-^hu ^> hcu in our wide trap limit (which has 
N — > oo) justifies the previous assertion that, at the phase transition, most particles in the 
cloud can be treated classically 

Now let's analyze the size of the region in which physics is non-perturbative at the 
transition. In our review of the uniform gas, we saw that physics becomes completely 
perturbative when p, c — fi ^> h a 2 /mA 4 . In an arbitrarily wide trap, the effective value 
of fi is p, — ^mu 2 x 2 . The condition for the existence of non-perturbative physics at the 
transition is then ^muj 2 x 2 ^ h 2 a 2 /m\ 4: , and the width of the non-perturbative region is 
R np ~ ha/mcuX 2 ~ A^ 1//3 a. Note that, even within this "non-perturbative" region, fluctua- 
tions with small wavelengths (^C are still perturbative. 

The relative volume of the non-perturbative region to the volume of the entire gas cloud is 
{Rnp /-Rcioud) 3 ~ ( a /0 3 - This means that non-perturbative contributions to the relation iV = 
/ d 3 xn(T, pL— V(x)) between N, T, and \x are suppressed by more than three powers of (a/l) 3 . 
It's more than three powers because, even in the relatively small non-perturbative regime, 
the dominant contribution to the density comes from typical particles, whose wavelengths 
are of order the thermal wavelength A ^> 1 /u and which can be treated perturbatively. This 
makes the total suppression (a//) 4 . In any case, the conclusion is that there is no obstacle 
at second order in a/ 1 to using perturbation theory to derive the relationship between N, 
T, and \x. 



III. n(T, /x) FOR A UNIFORM GAS AND ITS APPLICATION 

The second-order perturbative result for n(T, fi) can be easily extracted from an old 
second-order result of Huang, Yang, and Luttinger |T5|JT6|1 for the pressure of a uniform hard 
sphere gas: 

P = ^l [ U 5/2 (z)- 2 ^{U 3/2 (z)] 2 

/ n \2 ( oooooo 7 i+j+k \ \ 

+8 (a) ([^M]^M + ggg (i + t)(j + t)(W1/2 ) + -}. (3-D 

where 

z = (3.2) 

is the fugacity. Li n is the polylogarithm function, which for our purposes can be considered 
as defined in terms of its series representation, 

oo s 

Lin(*) = £-. (3.3) 

s=l 6 
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(a) (b) (c) 



FIG. 1. Perturbative Feynman diagrams contributing to the pre ssure P at (a) first order and (b,c) 
second order. Diagram (b) corresponds to the first (a/A) 2 term in (3.1) and diagram (c) to the second. 



We have independently rederived and verified this result for the pressure. For the sake 
of any readers who might find a derivation in the language of field theory [based on the 
Lagrangian ( |2.2| )1 a useful supplement to the original, we give the derivation in Appendix 
|A[ In the language of the imaginary-time field theory (|2.3|) , the perturbative Feynman 
diagrams which correspond to the first- and second-order terms in the pressure (|3.1|) are 
shown in Fig. |IJ. 

We now obtain n as dPj d\i: 



n 



1 f z l 2a (i + j)z i+ i 
A 3 I V i 3/2 A ^ (ij) 3 / 2 



+8 T + 



i+j+k 



+ 



(y)3/2 fc l/2 ({ + k)(j + k^ijk) 1 / 2 



+ 



(3.4) 



Here and henceforth, indices of sums (i,j,k) implicitly run from 1 to infinity. (Most of 
the terms above could be written in terms of polylogarithms, but the form shown is more 
convenient for the next step.) We emphasize that this is a perturbative expansion and is 
valid only in contexts where perturbation theory is applicable.^ In field theory language, 
the above result for n corresponds to the diagrams of Fig. 0. 

Now, use [L = p — V(x) and integrate over x, as in (|1.5|) . For the harmonic potential 



LTD , the integrals are simple Gaussian integrals, giving 



N 



ho . 



+1 



51 ,-3 



2a 



E 



§(^+J + ^) 1/2 



(v) 8/a (i + i) 1/a 
1 



+ 



1 



(y)3/2fcl/2 (i + k)(j + k)(ijk)V 2 



+ 



(3.5) 



2 The same expansion was incorrectly applied in Ref. jTJj to the problem of the first-order correction to T c 
for a uniform gas — a problem where perturbation theory breaks down. 
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(a) (b,) (b 2 ) 



(c) 



FIG. 2. Perturbative Feynman diagrams contributing to the number density n — (ip*ip) at (a) first order 
and (bi,b2,c) second order. The cross corresponds to an insertion of the operator ip*ip, whose expectation 
is being taken. 



where z = e"^ and 

is the geometric mean of the trap frequencies. 
In Section IV, we will discuss the expansion 



(3.6) 



/'c 



_(-\\0i _ 
He X +/ic 



,l( 2 ) ( - 



+ 



knT 



(3.7) 



of the critical value of /2 in powers of a. The ideal gas result is /2 C = 0. Recall from previous 
discussion that the second order term Ji^f 1 is not perturbatively calculable]^ We will later 
see that the second-order coefficient ft^ contains a logarithm, 



/4 2) 



A\n- + B, 
A 



(3.8) 



and that the coefficient A of the logarithm is perturbatively calculable, while B is not. 

In any case, we would like to insert the expansion ( |3.7| ) for Jjl c into the expansion ( |3.5| ) 
for N. The problem is a little more complicated than simply Taylor series expanding the 



3 For those readers who like to think in terms of Feynman diagrams, there is a hand-waving, heuristic 
argument to see why there is a non-perturbative 0{a 2 ) contribution to /2 C . The first-order contribution, 
given by the first diagram of Fig. |3|, contributes 0(a) to /i. Since this diagram is momentum independent, it 
can be absorbed into a renormalization of fi, so that one need not consider higher-order diagrams that contain 
it as a sub-diagram (a point relevant to the discussion of divergences at higher orders). At the transition, the 
second diagram of Fig. |^ contributes 0(a 2 ) times a logarithmic IR divergence, arising from the contribution 
where both loop frequencies are zero and the two loop momenta simultaneously approach zero. Three- 
loop contributions turn out to produce 0(a 3 ) times linear IR divergences, four loops produce quadratic 
divergences, and so forth. Suppose we heuristically cut off these infrared divergences at a momentum scale 
Air. The perturbation series then turns out to look like f3p, c = 0(a/X) + 0(a 2 /X 2 ) + 0(a 3 /A 4 Ain) + 

0(a 4 /A 6 Ajp) H , where we have ignored logarithms such as ln(AAiR). Imagine starting Air at some high 

value and then lowering it. The usefulness of the perturbative expansion will break down once we get to 
Air ~ a I A 2 , which is just the non-perturbative scale u discussed earlier. For this Air, all the terms in 
the series after the first become the same order, which is 0(a 2 /X 2 ). This suggests that this is the order 
of a non-perturbative contribution to /3/2 c . For a clean argument, however, one should instead refer to the 
analysis in the text. 
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individual terms of the sums of ( |3.5|) in £L, because such a procedure would lead to unregu- 
lated infrared logarithmic divergences at second order. We derive the small p, expansion of 
( P-5| ) in detail in Appendix IB], with the result 



N 



' k B r 

hu ho 
+ 



C(3) + 



A 



C(2)/i (1) -2^ 



-J p/2j3/2^ + jy/2 



|[^)] 2 + C(2)A (2) 



2/^ (1) E 



[i+j 



I 1 / 2 _ ;l/2 



+8E 



il/2 

-^-4C(|)/i (1) 



{ijfi 2 k 1 ' 2 + j + ky/ 2 (i + k){j + k){i + j + ky/ 2 k i / 2 i 



1 /a 

2 U 



(i) 



4C( 



In 



+ 



A 



(3.9) 



The logarithmic term at the end is the manifestation of the infrared logarithm just 
mentioned. In fact, at the critical point, the coefficient of this logarithm vanishes because 



/A 



4C0 



(3.10) 



A diagrammatic interpretation of why the logarithm vanishes is given at the end of Appendix 
|BJ. The first-order result (|3.10| ) for p^' can be derived using mean field theory, and a 
discussion in the context of trapped Bose gases may be found in the original first-order 
derivation of T c 0. We will rederive it in the next section, along with the second-order 
coefficient pf\ For the moment, though, let's use the known first-order result (|3.10| ) to 
solve for T c in terms of Inverting ( |3.9| ) gives 



T 



1 + Cl — + C 2 
A 



A 



+ 



A f 



(3.11a) 



ci 



c 2 



3C(3) 
C 2 



E 



C(2) 
3C(3) 



- 2C(2)C( 



-3.426 032, 



u (2) 

H'c ' 



(3.11b) 
(3.11c) 



C 2 = \c\ + 



3C(3) 



c(|) 2 + 2C(|)E 



'J 



-2E 



(, +J )l/2_,l/2_ J l/2 
^3/2 j 3/2 



+ 



>J 



(y)3/2jfel/2 y (i + j + fc)Va (i + fc)(j + jfe)(i + J + jfe)l/2 jfel/2 



21.4, 



(3.12) 



where 



' ^ \ V3 ^ho 
,C(3)J 



(3.13) 
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is the ideal gas result and 



An 



2nh 2 
\ mk B T 



2,71 



N 
X(3), 



-1/6 



^ho 



(3.14) 



is the corresponding thermal wavelength. The first-order result is the same as that found in 
Ref . @ .0 Results for the individual sums appearing above are listed in Appendix |C|. 



IV. /2 C (T) FOR A UNIFORM GAS 
A. Overview 

We'll now address how to relate the chemical potential fi appearing in the original 3+1 
dimensional theory ( |2.2| ) to the parameter r of the effective three-dimensional theory ( f2.5| ). 
The critical value r c of r can be extracted from lattice simulations of the latter theory 
which will then allow us to determine the critical value fi c {T) of fi. 

Effective theories, such as the three-dimensional 0(2) model, have long been used to 
describe long-distance physics at second-order phase transitions. Such use of effective theo- 
ries is often restricted to studies of universal quantities, such as critical exponents, because 
the relationship between the parameters of the effective theory and a more fundamental 
description of the system cannot be computed systematically. The situation is quite differ- 
ent for dilute Bose gases near the phase transition: the short distance scale A at which the 



long-distance three-dimensional effective theory description (2J5) breaks down is a scale at 
which the physics is perturbative (since A <C 1/w). One may therefore perform a perturba- 
tive calculation to relate r to /i, even though the long-distance physics at the transition is 
non-perturbative. 

Such perturbative matching of the parameters of effective theories with underlying 
physics has a long history in field theory. It has been applied in a number of problems, 



including lattice field theory |T9[, Bose condensates at zero temperature [EDI, relativistic 



corrections to non- relativistic QED pl~] , heavy quark physics p2| , ultrarelativistic plasmas 



2B| , and non- relativistic plasma physics |23| . For a general discussion, see also Ref. . The 



basic idea is to formally compute, in perturbation theory, some number of infrared physical 
quantities in both the effective theory and the more fundamental theory. By equating the 
results from the two theories, one can then solve for the parameters of the effective theory 
(to the order desired). 

The perturbative computations are performed using any convenient infrared regulator 
(though it must be the same regulator in both theories). The perturbation series for various 



4 The sum in ( |3.11b| ) is expressed in a slightly different form than in Ref. ||. The relation is that 



J \- (i + j) ,\- 1 

(i + i)V3 ■ 

where we've used i <-> j in the last step. 



j3/2^3/2^ _|_ j)l/2 2-j j3/2j3/2^ _|_ j)3/2 2-j j3/2jl/2^ _(_ j)3/2 ' 
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3+1 dim. theory of ip 


3 dim. theory of <p 




1 


1 


— > — 


ik + cu k 




> m > 


+fi 


—r 




8iia 


2u 




m 





TABLE II. Feynman rules, app ropriate for a matching calculation in a uniform Bo se gas, for the 
original 3+1 dimensional theory ( |2 3| ) of ip and the effective three-dimensional theory (2.5) of 4>. We have 
set % = fcs = 1. The variable fco represent the Matsubara frequency of the field, while ojk = k 2 /2m. At 
finite temperature, loop frequencies Iq are summed over the discrete values Iq = 2imT with n any integer. In 
dimensional regularization with the MS renormalization scheme, a factor of M e = (e 7B / 2 M/V^7r) c should 
also be associated with each 4-point vertex but has not been explicitly shown above. 



physical quantities will be badly behaved if one removes the infrared regulator since, in our 
case at least, the infrared physics is non-perturbative. But this bad infrared behavior will 
cancel out in the perturbation series derived for the parameters of the effective theory, and 
so one may safely remove the infrared regulator at the end of the matching calculation. This 
is a reflection of the fact that the difference between the effective theory and underlying 
theory has to do with short-distance physics, and short- distance physics is perturbative (in 
the cases where perturbative matching is applicable). 

The relevant distance scale of physics for the matching calculation is the short- distance 
scale A where the three-dimensional 0(2) theory breaks down. The corresponding energy 
scale is therefore of order h 2 /m\ 2 ~ k^T, which is simply the typical energy of particles in 
the gas. This scale is large compared to the size of the chemical potential at the transition 
(p.7|), which is of order (a/ \)k-sT. Therefore, for the purpose of doing a matching calculation, 
the chemical potential \i may be treated as a perturbation. In combination with the use of 
dimensional regularization, this turns out to be very convenient computationally. 

With /i treated perturbatively, the imaginary-time Feynman rules for the original 3+1 
dimensional action ( |2.3| ) are shown, for reference, in Table [TT| The analogous rules for the 
three-dimensional 0(2) effective theory are also shown. When discussing the evaluation of 
Feynman diagrams, we will always set H — &b = 1 in order to avoid cluttering up equations 
and discussions of conventions. We've specialized to the case of a uniform gas by taking 
V(x) = 0, and we've introduced the shorthand 



2m 



(4.1) 



We will use the notation k , l , p , ... to designate the Matsubara (imaginary time) frequen- 
cies associated with propagators with momenta k, 1, p, .... 

In our case, the short- distance length scale A -1 at which the three-dimensional theory 
breaks down is of order A, as we've discussed before. In principal, a long-distance 

effective theory can correctly describe physics at an infrared wavelength scale k -C A to any 
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desired order in k/A. However, as one pushes the description to higher and higher powers 
of k/A, one must add more and more corrections to the action of the effective theory, in 
the form of interactions that are more and more infrared irrelevant (in the sense of the 
renormalization group) — that is, interactions with higher scaling dimension. In our case, 
the long-distance physics scale of interest is the non-perturbative scale 1/u, and powers of 
k/A translate into powers of our expansion parameter uX ~ a/A <C 1. We shall discuss later 
why including such corrections, such as a |0*V0| 2 terms in the effective Lagrangian, would 
in particular not affect T c at second order. [We will also give a similar discussion of (0*0) 3 , 
which is a marginal operator in three dimensions.] For now, though, we shall simply ignore 
the issue and push ahead with the matching calculation. 

The action for a given effective theory can be written in a variety of equivalent ways by 
making field redefinitions, such as — > c0 or — > (l + eV 2 + - • ■)(/), etc. Our convention shall 
be to insist that the fields of the three-dimensional and 3+1 dimensional theories match up, 
to whatever order in k/A we are working, as 



^(0,k)= *|%(k). (4.2) 

This was our definition of in the more cavalier discussion in the introduction. The fre- 
quency ko of ip(ko, k) denotes imaginary-time frequency. So one of our matching conditions 
will be that the inverse Green functions match up as 



^(0,k) = S^ (4.3) 



to the relevant order in k/A. In the presence of interaction, this definition of might fix 
the normalization of the (V0)*(V0) term in the action (|2.5|) to be different from 1. Our 
three-dimensional effective theory should therefore be written in the somewhat more general 
form 



S = I d 3 ~ c x 



Z^(V0)* ■ (V0) + r e fj0*0 H — ^-(0*0) 2 + (higher-dimensional operators) 

(4.4) 

where can deviate from one beyond leading order. In principle, we need to determine the 
parameters r e ff, and w e ff (and any higher- dimensional operators, if they were required at 
a desired order) by matching. 



B. UV regularization 

Before starting on matching, we must first unambiguously define the parameters of our 
theories. The three-dimensional long-distance 0(2) theory (|2.5| ) is super-renormalizeable, 
but there are UV divergences associated with the 0*0 interaction. Diagrammatically, these 
divergences are associated with the graphs of Fig. || In order to give the coefficient r a well- 
defined meaning, we need to specify a regularization and renormalization scheme. By far the 
most convenient regularization scheme for perturbative matching calculations is dimensional 
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FIG. 3. Fundamental UV-divergent diagrams in the three-dimensional theory (2.5). Dimensional 
regularization automatically removes the linear divergence of the first diagram and regulates the logarithmic 
divergence of the second diagram as 1/e. 



regularization. We shall replace the number d = 3 of spatial dimensions by d = 3 — e, taking 
e — > at the end of the day. 

To define a finite, renormalized value of r, we will use the modified minimal subtraction 
(MS) scheme with a renormalization scale M. The theory is then 



S = I d 3 - e x 
with the relation 



^(V0)* ■ (V0) + r bar( #0 + M e -f 

o 



(4.5) 



1 fu\ 2 

between the bare coupling r bare and the renormalized coupling r m (M), and where 

p7e/2 

M = —= M. (4.7) 



[The factor of e^^/V^ in ( |4.7| ) is what distinguishes modified minimal subtraction (MS) 
from unmodified minimal subtraction (MS); the difference between the two schemes amounts 
to nothing more than a multiplicative redefinition of the renormalization scale.] 

The original 3+1 dimensional effective theory (|2.2|) is not renormalizeable and also re- 
quires UV regularization, and we will again use dimensional regularization. At second order 
in the interaction strength (the order relevant to our calculation), the only UV divergence 
is a well-known linear divergence associated with the second diagram of Fig. ^, which can 
be absorbed into a redefinition of the coefficient of the (ip*ip) 2 interaction. To relate this 
coefficient to the physical scattering length a, one needs to regularize the theory and then 
compute the zero-energy limit <r(0) of the 2-particle cross-section (at zero temperature and 
density), since a is defined by 87ra 2 = a(0) for identical particles. In dimensional regulariza- 
tion, however, the loop integral for the second diagram in Fig. [| vanishes at zero temperature 
and density, and so there is no second-order correction to cr(0). The coefficient of the quartic 
interaction therefore remains its tree-level value 2%h 2 a/m, as in (|2.2| ) frTfl . 



This property of dimensional regularization is a simple consequence of dimensional analy- 
sis. At zero energy (i.e. zero external momenta), the second diagram in Fig. ^ is proportional 
to the loop integral 

dL d A -n ^ 



(il + loi)(-Uq + ui) 
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FIG. 4. The 2-particle scattering amplitude at second-order in the interaction strength. (Note that 
there is no t-channel analog of the second diagram, since such a diagram vanishes in a non-relativistic 
theory.) 



All of the parameter dependence of this integral can be factored out by rescaling l Q by a 
factor of 2m. The rescaled integral 



dL d 3 ~H 



(ti + P)(-il + l 2 ) 



(4.9) 



has dimensions of (length) ~ 1+<E but no dimensional parameters to make up that dimension. 
The only dimensionally consistent answer is zero. In most regularization schemes other than 
dimensional regularization, there are still dimensionful parameters in the integral associated 
with cut-off scales, and the integral would not be zero. For example, if we regulated with a 
UV cut-off A on I, the integral would give a non-zero result proportional to A in d = 3. 



C. Matching of Z<p 



We want to calculate the critical value /i c to next-to-leading order [i.e. fi^ in the expan- 



sion Q3.7Q ]. One might expect this to require knowing the parameters r B g, and u e s of 
the effective theory to next-to-leading order. In fact, as we shall see, dimensional regular- 
ization organizes the calculation in such a way that we only need to compute r e g, which is 
three-dimensional analog of //. But let us briefly discuss the matching of anyway, as a 
simple warm up. 

The matching of Z^ is trivial because the first-order contribution to the inverse propa- 
gator, given by the first diagram of Fig. does not have any momentum dependence. That 



is, equating the inverse propagators of the two theories as in (|4.3|) gives 
k 2 1 

— + (k independent) + 0(a 2 ) = — Z^k 2 + (k independent) + 0(a 2 ) , (4.10) 
2m 2m 1 J 

where 0(a 2 ) indicates corrections that are formally second order in perturbation theory. So 

Z^ = l + 0(a 2 ). (4.11) 

In this paper, we will write 0(...) when displaying the full parameter dependence of a 
correction (except possibly for logarithmic factors) and write 0(...) when just showing the 
dependence on a particular parameter. So 32a 2 /X 2 = 0(a 2 /X 2 ) = 0(a 2 ). In matching 
calculations, where we are formally doing perturbation theory with IR regularization, 0(a n ) 
will just mean ra-th order in perturbation theory. 



15 



<7 





I I \ 






— • > • — 





FIG. 5. An two-loop example of a diagram contributing to — n^O) for a matching calculation, r is 
treated perturbatively. 

D. Matching of r 

To match r, take the k = case of matching (|4.3| ) the inverse Green functions: 



- 11 + n^(0, 0) = ^- [r bare + n^(0)] , (4.12) 

where II is the proper self-energy. 

We will use dimensional regularization to regulate the infrared divergences of perturba- 
tion theory, as well as the UV divergences already discussed. A well-known advantage of such 
use of dimensional regularization for matching calculations is that then every loop diagram 
contributing to n^(0) vanishes by dimensional analysis arguments similar to the one given 
in Section |1V B| . Consider, for example, the contribution of Fig. |5|. In the three-dimensional 
effective theory, this diagram is proportional to the loop integrals 

g (4.13) 



/V|i + q| 2 

which must vanish because there are no dimensionful parameters to make up the dimensions 
of the result. It is crucial here that there are no external momenta, that r may be treated 
as a perturbation for the purpose of matching calculations, and that loop integrals are never 
dimensionless in dimensional regularization. The upshot is that the matching condition 
( |4.12j ) becomes simply 

r bare = 2m[-/i + n^(0,0)]. (4.14) 

Another convenience of the vanishing, in dimensional regularization, of loop diagrams in 
the three-dimensional theory is that we need not keep track of the matching of u and if 
we're only interested in the matching of r. For instance, the one-loop contribution of the 
first diagram in Fig. |3] gives a contribution proportional to u in three dimensions, and so a 
second-order calculation of Ii^, would require a second-order determination of u, if it weren't 
for the fact that this diagram vanishes. 

Returning to the 3+1 dimensional theory, the diagrams which contribute to IT^ up to 
second order are shown in Fig. |], and all diagrams are to be evaluated at finite temperature. 
Diagram (a) gives the first-order contribution to IL^. It gives 
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FIG. 6. Diagrams contributing to —11^. /i has been treated perturbatively. 

<> = ^^-, (4.15) 
v m Tp tp + uj p 

where we introduce the shorthand notations 

t^l^M-JJ^, (4.16) 



Po *• PO 



where d=3— e is the number of spatial dimensions. The po sum can be performed by standard 
contour tricks,^ yielding 



n (a) _ »™ 

v m 

where n(u) is the Bose distribution function 



n(u p ) + l\ (4.17) 



»( w ) = iS^l- ( 4 - 18 ) 

The integral of a constant vanishes in dimensional regularization (again by dimensional 
analysis), and the integral of n(ujp) can be carried out in three dimensions to yield 



< a) ~C(§)+0(e). (4.19) 



4aT 
T 

At the critical point, diagrams (b) and (c) cancel each other and so need not be computed 
individually. The cancellation arises because the inverse susceptibility x^ 1 = — A 4 + n^(0, 0) 
will vanish at a second order phase transition. This condition is shown diagrammatically in 
Fig. ^ at first order in perturbation theory. As we've discussed in Section |11 C| , perturbation 
theory breaks down in the calculation of \x at second order, but the first-order relation of 
Fig. is therefore reliable.^ This relation implies that diagrams (b) and (c) cancel at second 
order: 



5 For example, see Section 25 of Ref. [p6| . 

You may wonder why we've discussed the reliability of perturbation theory here when we've already 
asserted that perturbation theory is valid for the matching calculation. The reason is that we're jumping 
ahead a little in order to streamline the calculation. The matching calculation can be done perturbatively 
because it involves only physics at the perturbative scale A, but the subsequent solution for r c (and therefore 
p, c ) cannot, since it involves physics at the non-perturbative scale 1/u. Since we are using a result about fi c 
to simplify our matching calculation at /2 C , we need to be careful. 
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= x' 1 = + LJ + Oia 2 ) 



FIG. 7. The vanishing of the inverse susceptibility at the phase transition, expressed in terms of 
diagrams at first order in a. (Unlike Fig. |l4], these diagrams are evaluated at the center of the trap and not 
integrated over x.) 

4 b) + = 0(a 3 ) at/i = /v (4.20) 



The sunset diagram 
We now turn to diagram (d), the sunset diagram: 

nl d) (0) = -- (— Yi t - • (4.21) 

^ 2\mJ TqTk (iq + uj g )(ik + uj k )[i(q + k ) +o; q+k ] 

We review in Appendix [B| how the loop frequency sums can be done with standard contour 
tricks, with the result 



^ 2 \ m J iqkl UJi — UJ q — LOj, 

(4.22) 



The symbol P.P. indicates the principal value prescription 

P.P. - = Re ] —, (4.23) 

where ± is an infinitesimal. This prescription removes the spurious divergence associated 
with uji — u q — Uk —>■ (but not also I ~ q ~ k —>■ 0), which is an artifact of this form of 
writing lT d )(0). (See Appendix 0.) We won't bother to explicitly write the P.P. in what 
follows. We note, as a side remark, that, in the language of time-ordered perturbation theory 
(in real time), the three terms of (|4.22|) correspond to the three diagrams of Fig. 
It's easy to see that the first term of ( |4.22[ ) vanishes, because it is proportional to 

n(u q ) n(uj k ) r n(u q )n(u k ) r n(u q ) n(u k ) 

^ . I„ , 1,12 „2 v> ~ . o„ 1, ' i 4 - 24 ) 



qk C<J q+ k — UJ q — UJ k -/qk |q + k| 2 — q 2 — k 2 Jqk 2q • k 

which vanishes by q — > — q (for any reasonable choice of regularization scheme). In Appendix 
0, we show that the last term of ( 4.22 ) vanishes as e — > in dimensional regularization. 



Only the second term of ([4.22 ) remains, giving 



(d)/M , fSira\ 2 f n(uj k )n(uji) 



n$»= — / V K/ V " +0(e). (4.25) 
v \ m J Jki uji - Uu-i - oo k 
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FIG. 8. Time ordered perturbation theory diagrams corresponding to the three terms of (4.22). Time 
flows from left to right, and the dashed line represents the zero external momentum (and zero frequency) in 
diagram (d) for —11.0(0,0). Each incoming momentum (other than the zero one) is associated with a Bosc 
distribution factor n. 



Subtracting divergences of the sunset diagram 



Because n(uj p ) — > (j3uj p )~ l oc p~ 2 as p — > 0, the above integral ( |4.25 ) has a logarithmic 
infrared divergence associated with k ~ I — > 0. Unfortunately, the full integral in ([4.25 ) is too 
complicated for us to do in arbitrary dimensions, which would be the most straightforward 
way to apply our chosen regularization scheme, dimensional regularization. It's convenient 
to instead explicitly isolate the divergent IR behavior by rewriting 



ni d) (o) 



87ra\ 2 
m J 



kl 



n(u k ) n(ui) - n (oj k ) n {ui) 



n (uj k ) n (ui) 

kl UJl - LU k -\ - UJ k 



+ 0(e) 



where 



nolo;) 



1 



(4.26) 



(4.27) 



The second integral vanishes in dimensional regularization for the usual reason: it is pro- 
portional to 



1 



/ki k 2 l 2 {P - |k- 1| 2 - k 2 ) ' 
which contains no dimensionful parameter to make up its dimensions. So 
r(d)/n\ f8ira\ 2 f n(u k )n(ui) - n (uj k )n (ui) 



(4.28) 



C ; (0) 



m 



ki 



— ^k-l _ 



0(e). 



(4.29) 



The above integral is infrared convergent and, if it weren't for the fact that we've now 
introduced a UV divergence associated with k ~ I — > oo, we would be able to set d = 3 in 
that integral and ignore regularization issues. 

To continue, it is useful to understand another way to interpret the infrared behavior 
represented by the last term of ( 4.26 ): 



/87ra\ 



V m 



n (u k )n (uji) 
ki uji - u; k _i - uj k ' 



(4.30) 
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As discussed before, infrared physics is dominated, in imaginary time, by the zero-frequency 
mode of the field if). Iq turns out to be the q = k = piece of the original frequency sums 
( |4.21| ) representing diagram (d). A quick way to see this is to note that the integrand in I 



above is the high temperature limit {(3 — > 0) of the original integrand in ( 4.25Q . But, if one 



goes all the way back to the original imaginary-time frequency sums ( [4.21 ), the integrand 
there is proportional to 

w (*9o + u q ) {ik + u; k ) [i(q + k ) + cj q+k ] ' 



with go and ko of the form 2irnT. Only the go = — piece survives in the (3 — > limit of 
this integrand, and this establishes the correspondence. 

It is important to note that non-zero frequency modes do contribute to diagram (d) 
even in the infinitely high temperature limit, because the limit does not commute with the 
integration over spatial momenta q and k. However, in our analysis so far, we have not yet 
performed the q and k integrations, and it is okay to take limits of integrands to see the 
correspondence of I with the q = k = piece of diagram (d). 

The upshot is that the infrared piece I that we isolated from diagram (d) is proportional 
to the same diagram evaluated in a purely three-dimensional theory: 

T 2 /8vra\ 2 f 1 A ~ , /8vra\ 2 r 1 /in . 

— 4m T ( — — ) /-___. (4.32) 



2 V m J Jqk a;gC<j fc u; q+k V m / J q kg 2 /c 2 |q + k 

This diagram is logarithmically divergent in both the infrared and ultraviolet, just as the 
original expression ( |4.30| ) for Jo, and it vanishes in dimensional regularization. The UV 
divergence of our current expression ( 4.29|) for diagram (d) came from the UV divergence 



of Iq. To isolate this UV divergence, we'd like to isolate a term that (i) has the same UV 
divergence as Iq, (ii) is analytically computable in dimensional regularization, and (iii) is 
infrared convergent (since otherwise we'll just re-introduce an infrared divergence when we 
isolate it). Something which satisfies all these requirements is the same integrals (|4.32|) of a 
three-dimensional theory as above but with mass terms to cut off the infrared: 

T 2 /8vra\ 2 f 1 

In = 



2 V m / 7qk (u q + AT) {u k + AO (cj q+k + A/")] 

f^V / ± , _ , (4.33) 



m J V qk (g 2 + A4 2 )(F + A4 2 )(|q + k| 2 + A4 2 ) 

where M = Ai 2 /2m is an arbitrary frequency scale. Our strategy will then be to rewrite 
our current expression ( |4.29| ) as 



7ray r n(u k ) n(ui) - n (uj k ) n (uJi) ^ 

Jk\ LOi — Wt-i — UJh 



n»(o) . 

m / Jki uji- cu k _i - u k 



I N + 0{e). (4.34) 



To put the first 1^ term in a form similar to the integral shown explicitly in (|4.34[) , one 



may replace u by uo + J\f in our early discussion of Ili^ (0) and take the (3 — > limit in all 



integrands to get the following analogy to (|4.22j ): 
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IN 



8na\ 2 f [n (uj q + Af)n (u k + Af) - 2n (uJ k + Af)n (ui + M)\ 



■ 



qkl 



We then obtain 



ni d) (o) 




(w, + M) - (u q + A/") - (u k + Af) 

x(27r) d 5 (d) (l-q-k) + C(e). 

n(a? fc ) n(ui) - n (uj k ) n (uji) n (uj k + Af) n (uji + Af) 



(4.35) 



+ 

Wo(^-i + A/") n (uj k + A/") 



ki 



AA 



- /at + O(e), 



AT 



(4.36) 



The first integral (with its implicit P.P. prescription) is now both infrared and ultraviolet 
convergent and can now be evaluated in exactly d=3 dimensions. The second integral is 
convergent as well. So fix d=3 in these integrals, scale out the parameters, and do the angular 
integrations using the P.P. prescription. This puts the integrals in a form appropriate for 
straightforward numerical evaluation. The result for the dimensionally regulated integral 
O) for I N is 

/ 

Jal 



(An) 



1 , M 
h In— — 

2e 3M 



qk (q 2 + M 2 )(k 2 + M 2 )(\q + k| 2 + M 2 ) 

which can be extracted from the general d result of Ref. |27| or the e 
28fl . Putting everything together, 

327ia 2 T 



1 

+ 2 



<>(<>) 

where C% is the numerical constant 
1 



2e 



ln(367rA^) 2 



tt* Jo 



dk dl 




+ ln(MA) + d 



kl 

ol 2 



+ 0(e), (4.37) 
analysis in Ref. 

(4.38) 



kl 



J In 


fc-Z 







(F+A r )(/ 2 +A r ) 



In 



Af + 2k(k-l) 



Af + 2k{k + l) 



1 



71"^ JO 



<ig eZ& 



(g 2 + AT)(A; 2 +AT) 



In 



Af + 2qk 



Af-2qk 



which is independent of the choice of the dimensionless number Af = /3Af. 
evaluation of the integrals gives 



-0.54410. 



(4.39a) 
Numerical 

(4.39b) 



(Since completion of this work, a somewhat more compact formula for C\ has been derived 
in Ref. 0].) 

We should mention that it is possible, at a formal level, to turn the original unregulated 
integral of (|4.25|) into a double sum, similar to the sums appearing the earlier formula fl3.4| ) 
for the density, by using methods similar to those reviewed in Appendix However, the 
infrared divergence of 11^(0) would manifest as i,j—>-oo divergences of these sums. We 
found it easier to handle the infrared issues in the integral form than in the summation form. 
This is the only reason why our treatment of 11^(0) superficially looks so dissimilar, in final 
form, to our treatment of pressure and density in Section [Hj . 



21 



1. Final result for r 



Combining our results for the pieces of IL^O) with the matching formula (|4.14j ) for r, we 
obtain 



_ f 4a w ,. 327ra 2 
-2mfi + 2mT — C(|) + 



2e 



+ ln(MA) + d 



+ 0(a 3 ) + 0( f i- f i c ). (4.40) 



Comparing to the expression (|4.6| ) for the MS definition of r, and using the leading-order 
result (12.61) for u, 



Tms(M) = -2m// + 2mT I —((D 



32vra 2 



\n{MX) + C 1 +0(a 3 ) + 0(fi-ft c ). (4.41) 



E. Final result for /i c 

We can now solve Q4.41D for the coefficients in the expansion 



I'c 



r,w: 



,( 2 ) 



+ ■ 



of /i c . The first-order result, well known in the literature, is 



M 1) =4C(| 



5 ) 

2^' 



(4.42a) 



(4.42b) 



and is a simple consequence of the vanishing susceptibility as depicted in Fig. [7|. The second- 
order coefficient is 



p,® = 32tt 



ln(MA) + Ci - 72tt 



2 r Cj Ms(M, u) 



(4.42c) 



where r Cj Mg(M, u) is the critical value of r-^g for a given choice of coupling u and renormal- 
ization scale M. The only dimensionful scale of the three-dimensional theory at its critical 
point is u, and so one should pick the renormalization scale M of order u. Note that the 
critical value r c ^g is then proportional to u 2 by dimensional analysis. Taking M = u/3 for 
definiteness, and because that was the choice made in presenting lattice simulation results 
in ref. ITU, we have 



/4 2) = 32tt 

where the dimensionless constant 



ln (^) +Cl _ 72 ^ 



n = 



r c<m (M=u/3, u) 



(4.43) 



(4.44a) 



is non-perturbative and must be extracted from simulations of the three-dimensional effective 
theory (|2.5| ). The simulation result is, from Eq. (1.5) of Ref. @, 



K = 0.001920(2). 



(4.44b) 
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V. FINAL RESULT FOR T c 



We can now combine the second-order result ( J4.43| ) for pL c with our earlier expression 
for T c to obtain 



T C = T 



with 



ci ~ -3.426 032, 



(5.1a) 



(5.1b) 



32ttC(2) 
3C(3) 



(5.1c) 



32ttC(2) 



ln(327r 2 ) + Ci - 72tx 2 K 



-155.0 



(5.1d) 



and with T , A the ideal gas quantities given at the end of Section [TII|. The constants ci, C 2 , 
Ci, and 1Z are given by Eqs. (|3.11b| ), (|3.12|) , (|4.39| ), and ( f4.44|) respectively. All have been 
computed perturbatively, except for TZ, which is the non-perturbative information extracted 
from lattice simulations. 



VI. YET HIGHER ORDER CORRECTIONS 



We have based our discussion on the 3+1 dimensional theory ( |2.2| ) of ip and the effective 
3-dimensional theory ( |2.5|) of the zero-frequency Matsubara modes. Both of these theories 
are approximate and have corrections which we have ignored, claiming them to be higher 
order than the order of interest. In this section, we will briefly discuss the nature of those 
corrections. 

Let's begin with the original 3+1 dimensional theory ( j2.2| ) of ip. Among other things, this 
theory ignores (a) the energy-dependence of the low-energy atomic scattering cross-section, 
and (b) the effects of 3-body collisions. Braaten, Hammer, and Hermans [35] give a nice 
discussion of how to systematize the corrections to the low-energy 3+1 dimensional theory, 
discussing interactions that are progressively more and more irrelevant at low energies. The 
most important such corrections are to supplement the Lagrangian ( |2.2j) by the additional 
interactions 



ith 2 a 2 r s ,„ 



2m 



\V(ip*ip)\ 2 - B(ip*ipy 



(6.1) 



B parameterizes the amplitude for 3-body collisions. r s is the effective range of the 2-body 
scattering problem and parameterizes the linear term in the energy dependence of the cross- 
section at low energy. The importance of the r s term grows with energy, which turns out 
to mean that its leading effect on the critical temperature or the chemical potential is not 
infrared dominated and can be treated perturbatively. 
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The parametric size of the leading-order effects of these corrections can be estimated in 
a very simple way by comparing them to the usual quartic term (ip*il)) 2 . At leading order, 
the effects of (ip*ip) 2 on the quantities computed in this paper [/i c for a homogeneous gas 
and N(T C ) for a trapped gas] were dominated by momentum scales of order k ~ 1/A (as 
opposed to the infrared scale u). Relative to the a(ip*ip) 2 interaction, one would expect 
that the leading-order effects of the a 2 r s \ V(tp*ip)\ 2 interaction of (|6.1|) should therefore be 
suppressed by ar s k 2 ~ ar s /X 2 . Near the transition, this is down by two powers of the 
typical inter-particle separation I ~ n^ 1 ^ ~ N^ l ^a^ discussed in the introduction (since 
A ~ I at T c ) , whereas the second-order effects computed in this paper are down only by one 
power, compared to the leading-order effect of interactions. Again relative to the (ip*tp) 2 
interaction, one would expect that the 3-body (if)*if)) 3 interaction of (|6.1|) is down by a factor 
of (Bm/h 2 a)ip*ip ~ {Bm/Ti 2 a)n, which is down by three powers of / since n ~ / 3 . The 
moral is that corrections to the original 3+1 dimensional Lagrangian (|2.2|) do not matter 
for a second-order calculation of T c for a dilute trapped gas, that the result at third order 
would depend on the effective range r s and not just the scattering length a, and that the 
result at fourth order would depend on the 3-body scattering rate as well. 

One can also verify the above analysis by a consideration of the leading-order diagrams 
involving a given correction from ( |6.1| ). Fig. |S| shows diagrams contributing to /2 C and Fig. [II] 
those[] contributing to n(T,fi). As an example, the diagram of Fig. |9]a gives a contribution 
to the chemical potential proportional to 



(6-2) 



a 2 r s r k 2 
hm T ik + loi 

The diagram is not dominated by infrared momenta, and so the perturbative treatment is 
justified. The dominant wave numbers are k ~ 1/A, as claimed above, corresponding to 
energies k-^T and frequencies uj^ ~ k B T/h. The result is that 

a 2 r s k^Tk 5 a 2 r s 
5 I X ~ * vT fc BJ, (6.3) 

where the k 5 comes from the k 2 d 3 k in the integral. Compared to the leading-order result 
0(k B Ta/ A) for the chemical potential, (|6.3|) is down by 0(ar s /X 2 ), just as we argued more 
simply above. 

Finally, even ignoring corrections to the original 3+1 dimensional theory, there will still 
be corrections to the effective 3 dimensional theory ( |2.5| ) of the zero modes. One might 
worry in particular about a (0*0) 3 interaction between the zero modes, which is a marginal 
interaction in three dimensions. Such an effective interaction can be induced by diagrams 



such as Fig. |TT| in the 3+1 dimensional theory, where the external lines are zero-modes and 
the internal lines are non-zero modes. However, the non-zero modes are infrared insensitive 
and are dominated by frequencies of order k-^I '/h and momenta of order k ~ 1/A. Power 
counting Fig. then gives an interaction in the effective 3 dimensional theory of order 



7 The effect of r s represented by Fig. [Toj a has been considered historically in discussions of AT C for a 
homogeneous Bose gas ]3(|. Those discussions completely missed the dominant contributions to AT C . They 
also did not use the more general language of effective ranges but implicitly used Born approximation to 
express r s in terms of the 2-body potential. 
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(a) (b) 

FIG. 9. Leading-order corrections to the 11^,(0) (and hence the determination of p c ) due to (a) the 
effective range and (b) the 3-body scattering terms of ( |6.l| ) . The dotted line represents the momentum flow 
the \V(4>*ii))\ 2 vertex. 





(a) 



(b) 



FIG. 10. As Fig. pJ but showing corrections to the density n{T,jj). 
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FIG. 11. An effective (<p*4>) 3 interaction of zero modes generated by a loop of non-zero modes. 



5^ ~ w 3 A 3 (0*0) 3 , (6.4) 



where the w 3 can be understood as arising from the three vertices in Fig. |y] and then the 
A 3 from dimensional analysis based on the dominant momentum scale. 

Now consider the effect of the vertex (|6.4j) on the infrared physics at momentum scales 
p ~ 1/u, to which the 3-dimensional effective theory is intended to be applied. At that scale, 
the (0*0) 2 interaction can no longer be treated perturbatively and, by dimensional analysis, 
the fluctuations in are of order u 1 ^ 2 . The relative importance of the (0*0) 3 term at the 
infrared scale p ~ 1/u is then 

^A 3 (0*0) 3 7/ 3 A 3 « 3 ^ 

u(M)* ~ uX ^A 3 " (6 ' 5) 

The contributions of the effective (0*0) 3 operator is therefore down by three powers of 
/ _1 ~ A -1 compared to those contributions we have included in this paper. Other corrections 
to the three-dimensional theory are similarly suppressed. 



VII. HOW WIDE IS A WIDE TRAP? 

We have assumed throughout that the trap is arbitrarily wide (u x , u y , oj z — ► with 
Nuj x uo y uj z fixed). We will now take a moment to explain parametrically how wide "wide 
enough" is for our second-order results to be valid. Our second-order result for p, c depends 
on non-perturbative physics near the center of the trap, and we treated the trap as flat 
over the wavelength 1/u of such physics. The trap must therefore be wide enough that this 
wavelength fits comfortably inside the region of the trapped gas that is nearly critical, whose 



size we labeled R np in Section flD 



First consider the case u x ~ u y ~ u z . Using Table |, the condition 1/u <C -R n p can be 
translated into I A^ 1 / 6 a. Combining this with the basic diluteness assumption a I of 
our analysis, we then require 

a < / < N 1/6 a. (7.1) 

This shows only the parametric dependence, and we have made no attempt to estimate 
numerical factors. 

For a very anisotropic trap, the strongest constraint will come from requiring the nar- 
rowest direction of the near-critical region to be larger than 1/u. Let u; max be the largest of 
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I — I 

5T 

c 



FIG. 12. A schematic depiction of finite-size effects on the BEC phase transition in the ideal-gas limit. 
The dashed curve indicates the infinite-volume transition. 



y- 



and oj z . Repeating the analysis of Section IIP then gives the corresponding value of 



_R np in that direction as ha/muj mSLX X 2 ~ -/V^auv/^Wx- So the condition is 



< / < N 1/6 i 



0J n 



(7.2) 



These constraints may be translated into other variables using I ~ n -1 / 2 ~ A ~ iV~ 1//6 ah - 

This condition on the size of the trap can also be summarized as a comparison of the 
uncertainty in the value of T c due to finite size effects vs. the resolution with which we 
have computed T c in our second-order formula (|5.1| ). Finite size effects round off the non- 



analyticity of the infinite-size transition, as depicted in Fig. [12|. A standard result from the 
literature is that, below the rounded transition, finite-size effects create the appearance of a 
transition temperature shifted by M,WN 



* y c „ C(2) g 1/3 
T 2C(3) 2 / 3 ^ ho 



(7.3) 



from the infinite- volume value (in the ideal gas approximation), as depicted in the figure. 
Here, uo is the arithmetic mean 



UJ X + UJ V + UJ Z 



(7.4) 



The condition that this finite-size effect on the transition be small compared to the rela- 
tive 0(o 2 /Aq) correction to T that we have presented in (|5.1| ) happens to be the same, 
parametrically, as the right-hand condition in ( |7.2| ). 
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FIG. 13. An additional perturbative diagram contributing to the pressure at second order. The 
fat black dot represents the one-loop (ip*tp) 2 renormalization counter-term, which vanishes in dimensional 
regularization. 



VIII. CONCLUSION 

The relative size of the second-order effect in our final result ( |5.1| ) for T c obviously 
depends on the diluteness of the gas and the value of the scattering length, which will vary 
from experiment to experiment. However, just for fun, let us put numbers to the size of the 
effect for one specific experimental study of T c that has appeared in the literature. The 1996 
experiment of Ensher et al. [[32] found AT C /T = —0.06 ± 0.05 for dilute gases of roughly 
N = 40, 000 87 Rb atoms in the F = 2 hyperfine state, trapped with u z = 373 Hz, u z = 27rv. 



z; 



and Ul x UJy UJ Z /V8. The relevant scattering length is a = (103 ± 5) ao ||33|| , where 

ao = 0.0529177 nm is the Bohr radius. (See also Ref. ||34|| .) These parameters correspond to 
a/Ao — 0.016. For an arbitrarily wide trap, this would translate into a first-order correction 
to T c of roughly —5.4% and a second-order correction of roughly +0.9%. For the actual 
trap, however, the corrections ( |7.3|) due to finite-size effects are roughly —2.4%. The fact 
that this is larger in magnitude than the second-order correction leads us to suspect that 
this particular trap may not be wide enough for the second-order result to be trusted, as 
was discussed in Section IV11. 
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APPENDIX A: FIELD THEORY REDERIVATION OF P(T,fi) 

The two diagrams which contribute to the pressure at second order in perturbation theory 
were shown in Fig. |](b,c). Technically, there is also a third diagram, Fig. [H| which involves 
the one-loop (ip*ip) 2 counter-term (represented by the fat dot) for renormalizing the linear 
UV divergence of zero-temperature, zero-density 2 — > 2 scattering at second order. However, 
we shall use dimensional regularization, for which this counter-term vanishes, as discussed 
in Section [IV B| . The Feynman rules are given in Table |T|. They are the same as in Table |Tl 



except that we are not treating the chemical potential /i as a perturbation in this context. 
We work in units where h — &b = 1. 
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TABLE III. Feynman rules appropriate for standard perturbation theory in the 3+1 dimensional theory. 



1. The basketball diagram 

Let's start with diagram (c) of Fig. |I[ The corresponding contribution to the pressure 
P= (pv)- 1 \nZ if] 



P5 : 



Po+qo-ko-k) 



(27r) 3 5( 3 )(p + q- k - 1) 



m ) Tpqkl (ipo + ujp) (iq + ui q ) (ik + uJ k ) (Ho + &i) 
where we have introduced the shorthand notation 

k 2 

u k = U k - j2 = fl, 

2m 



(Al) 



(A2) 



and where the summation- integration sign is defined in ( |4.16|) . We now use standard tricks to 
evaluate the frequency sums.0 Specifically, rewrite the frequency Kronecker 6 as an integral 
of exponentials, and factorize the expression into independent sums: 



-ictpo 



-ictqo 



m 2 Jo Tp ip + uj p 7q iq + ui q jk ik + uj k % il + uji 



\-iako 



— (27r) 3 5 (3) (p + q-k-l) 



Then we use the standard frequency sums 



TJ2 



Pu 



tJ2 



e -iap 

ipo + w 

■,+iapo 



n(u) e c 



n[uo) e 



{/3-a)oj 



for < a < (3. The a integration is then trivial, yielding 

e /3(CJ p +Cd q ) _ e /3(ui h +uii) 



(A3) 

(A4a) 
(A4b) 



Pr 



^ 2 a 2 



m z Jpqkl 



n(u p ) n(u q ) n(Cb k ) n(ui) 



u p + u q -u k - UJi 



(27r) 3 (5 (3) (p + q-k-l). 



(A5) 



8 VP C corresponds to Q'/Q^ ofRef. ]l5[ . The contribution VPt, from diagram b corresponds to [Q" / Q^)- 

1[ S (D/ S (0)]2. 



See, for example, Section 5.5.1 of Ref. 
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Note that the zero of the denominator at Cj p +uj q = tik+uJi is canceled by a corresponding zero 
of the numerator. However, it will be useful to split the integral into pieces that individually 
lack this cancellation, and so it is useful to first introduce a redundant principal part (P.P.) 
prescription in (|A5|). Making use of the identity 

(A6) 



n(uo) e/** = n(ui) + 1, 
expanding terms, and permuting integration variables, we can rewrite 



as 



Pr 



32ttV 

m 2 
32vr 2 a 2 

m 2 



The term 

Pc,2 



P.P. 



n[uj q ) n{Uk) n{uji) 



\n(£j p ) n(u q ) 



pqkl 



P.P. 



pqkl 



32vr 2 a 2 



77T 



Up + UJq - UJ k - UJl 

n(u q ) n(u k ) n(ui) - \n(uj p ) n(uj q ) 

LU p + UJ q - LU k - UJl 



(27r) 3 <5 (3) (p + q-k-l). 
(27r) 3 (5 (3) (p + q-k-l). (A7) 



P.P. 

pqkl Ulp + UJq 



^n{Up)n{u q ) 



U) k - UJl 



(27r) 3 5( 3 )(p + q-k-l) 



(A8) 



involving just two n's, has a linear UV divergence associated with k, I — > oo with p and q 
fixed. This is the divergence that is canceled by the counter-term diagram of Fig. O for 
generic regularization schemes and which dimensional regularization will simply ignore. In 
fact, the entire term P c 2 simply vanishes in dimensional regularization, which can be seen by 
doing the k and 1 integrations explicitly in d spatial dimensions. Defining s = k — |(p + q), 



(2^g(p + q-k-l) 
P.P. = — m 

kl UJ p + UJ q — UJk — UJl 



d d s 
(2tt)' 



P.P. 



o 1 I I 

s 2 ~ 4IP - q| 



(A9) 



It's convenient to re-express the principal part in terms of infinitesimals, using (|4.23| ), before 
doing the s integration. The integral ( |A9| ) then yields 



mT (l - I) 



(47r) d / 2 



Re 



mT ( 1 — 



-ilp-q| 2 + ^o ± 

(d-2)/2 



(d-2)/2 



ilp - q| 



(Aix) d l 2 

Analytic continuation to d — 3 yields zero: 

(27r) d 5( fi )(p + q-k 



cos 



(d-2)7r N 



(A10) 



P.P. 



1) 



kl 



UJp + UJq - UJ k - UJl 



0(e), 



(All) 



where e = 3 — d. To conclude that the contribution ( |Aq ) to the pressure vanishes in 
dimensional regularization, one must also check that the final p and q integrals with the 
(9(e) integrand do not diverge, since divergences could possibly generate a 1/e singularity 
to cancel the 0(e) behavior of the integrand. However, the UV is cut off by the distribution 
functions n(uj p ) and n(uj q ) in ( |A8|) , and so this is not an issue. 
We are left with only the term of (|A7|) with three n's: 
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Pr 



32n 2 a 2 



P.P 



pqkl 



n{u q )n{u k )n{ui) 

UJp + UJ q ~UJ k - OJi 



(27r) 3 5 (3) (p + q-k-l) 



(A12) 



This reproduces eq. (A15) of Huang, Yang, and Luttinger flI5| . Since their subsequent dis- 
cussion of evaluating this integral is somewhat telegraphic, we will present our own method. 
First, expand the distribution functions n(Cj) = n(u — ji) in powers of fugacity z = exp(/?/i): 

327r 2 a 2 oo oo oo 



P 



EEE 

a=l b=l c=l 



a+b+c 



P.P. 



pqkl 



e -a/3q 2 /2m e -b(3k 2 /2m e -c/3l 2 /2m 



(27r) 3 (5 (3) (p + q-k-l) 



Rescaling all momenta by y/3/m to make them dimensionless gives 

8a 2 T 00 00 00 



A 5 



EEE- 

a=l fe=l c=l 



i-abci 



(A13) 



(A14) 



' abc 



(27T) 



9/2 



P.P. 
pqkl 



e -ag 2 /2 e -bA: 2 /2 e -c/ 2 /2 

i(p 2 + g 2 -P-/ 2 ) 



(27r) 3 (5 (3) (p + q-k-l). 



(A15) 



For the sake of justifying later manipulations, it is convenient to introduce a redundant 
exp(— + p 2 ) convergence factor into the integral defining I a b c . We will evaluate l a bc by 
representing the energy denominator and the 5 function as integrals of exponentials. Using 
the infinitesimal version ( f4.23| ) of the principal part prescription, we write 



Lbc = (27T) 9 / 2 Re / e -0V e -V/2 e -^/2 e -d 2 /2 
J pqkl 



r 1 00 r 
f d\ e -^W-k 2 -P^)X/2 J d 3 x e ix.(p + q-k-l)_ 



(A16) 



The p, q, k, and 1 integrations are now simple Gaussian integrals, yielding 



I abc = (27r)- 3 / 2 Re J dX e i0+x J d 3 x (0+ + A)- 3 / 2 (a + A)- 3 / 2 (6 - A)- 3 / 2 (c - A)" 3 / 2 

' x 2 f 1 1 1 1 \1 

x exp — H + H . (A17) 

I 2 V0+ + A a + A b-X c — A/J v ; 

The exp(zO + A) prescription is now redundant and can be dropped. It's also convenient to 
change integration variables from A to A+0 + in order to remove the remaining + prescription 
from the integrand (noting that a, b, c 7^ 0). The x integral is Gaussian and yields 



/■ioo+0+ 

I abc = Re / dX 



abc + 2bcX - (a + b + c)A s 



-3/2 



The final integral is straightforward and gives 

1 



' abc 



{a + b)(a + c)(abcy/ 2 ' 
The final result for this contribution to the pressure is 



y a+b+c 



Q n 2 r p 00 00 00 
Oil J. % ^ v -v v -v 

C = h ti k (a + b)(a + c)(abcy/ 2 " 



(A18) 



(A19) 



(A20) 
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2. The three-circle diagram 



The other second-order diagram, Fig. |l]b, is trivial in comparison. It's contribution to 
the pressure is 



Ph 



1 



Sua 



m J Tp (ip + (j p ) 2 \T iq + uj q 



One of the summation-integrals is 
1 



n[u q ) + 2 



m 



n(u 



3/2 



The other is easily obtained by differentiating with respect to /i: 



1*3/2(2) 



p (ip + dip) 2 



P 



ni 



3/2 



Lii/2^)- 



So, 



A 5 



-[Li 3 / 2 (2;)] Ui /2 (z) 



(A21) 



(A22) 



(A23) 



(A24) 



Putting P h and P c together gives the total second-order contribution to the pressure 
which appears in (|3.1|) . The first-order contribution of Fig. |l|a is easily evaluated in a similar 
manner. 



APPENDIX B: SMALL p, EXPANSION OF N 
1. The expansion 

Consider the (a/ A) 2 term in the expansion Q3.5p for N. First consider the term propor- 
tional to 

^i+j+k 

E (zj) 3/ 2A ,l/ 2(z+J + A;) l/2- ( B1 ) 

Because of the explicit factor of a 2 , one might naively think one could use the order a result 
z c± 1 for z. But this would give 

1 

§(u) 3/2 fc 1/2 (« + J' + fc) 1/2 ' (B2) 

which has logarithmic divergences associated with k — > 00 with 2 and j fixed. We can isolate 
these divergences by rewriting the original sum as 



^i+j+k 

g (zj) 3/ 2A; l/2 (z+J + fc) l/2 =E 



(i J f/ 2 k l / 2 {i + j + ky/ 2 (i]fi 2 k 



+ S«F^- (B3) 
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We can now safely set z to 1 in the first sum on the left-hand side. The second sum is easy, 

2' 



giving — £(|) 2 ln(l — z). The small // result is then 



^i+j+k 



E (ifi3/2 k l/2(i + j + k) l/2 



E 

ijk 



1 



{ijf^k 1 / 2 ^ + j + k) 1 / 2 {ijf' 2 k 



-C(|) 2 ln(-/3/i) + 0(/2). 



(B4) 



The sum associated with the a/ A term of ( |3~5| ) must be expanded to first order in /i, 
where it suffers a similar problem. Naively, 



gi+j 1 (, _l 7 'U/2 

ij (ij) 3/2 (i + j) 1/2 ^{ijfl'ii + j) 1 ' 2 {ijY' 2 



(B5) 



The second term has logarithmic divergences associated with (i) i — > oo with j fixed, and, 
symmetrically, (ii) j — > oo with i fixed. Proceeding as before, we can isolate the divergent 
behavior by writing 



E 



z i+j 



{ij) 3/2 {i+j) l/2 



E 



yi+j 



(ij) 3 / 2 (i+j) 1/2 i 2 P /2 i 3/2 j 2 



+ E 



ij 



+ 



^2^ 3/2 ^3/2^ 2 



(B6) 



In the first sum, we can now safely replace z by 1 + f3fi + 0(/z 2 ), and the second sum gives 



2C(|)Li 2 (z) = 2C(f ){C(2) + [- H-m + 1] pp) + 0(fi 2 ). 



(B7) 



The final result for the expansion is then 

1 

(ij)3/2(i +j )l/2 = N 



Et^ 



E 



;j ( , j) 3/ 2(z + j) T7I + 2 C(|)^[-ln(-/3/i) + l] 

'j + j)l/2 _ ,1/2 _ ,1/2 



(y)S/2 



+ 0(/i 2 



(B8) 



The last thing we need is the expansion of the sum in the order a term of (|1 . 5|) , which 
is just 



Li 3 (z) = C(3) + C(2)/3/i + | [- ln(-#Z) + §1 (/?/i) 2 + 0(^ 3 ). 



(B9) 



Combining the expansions ( |B4|) , ( p8|) and (p9|) of the sums with the expansion Q3.5|) of N, 
we obtain 
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N 



Hu- 



bs , 



+ 



C(3) + 



2a 



C(2)/3a*-xE 



A ^i3/2 j3 /2^ +J )l/2 



A U 



(i + j) 1 / 2 - i 1 / 2 - j 1 ' 2 4a 



1 



p/2j3/2 
1 



A 



ftiC(§) 



+8 LJ § '{ijfW 2 \{i + ] + kyi 2 1 (j + jfe)(j + ^(1+^ + ^)1/2 



-|(^-^C(D) m(-/3/2) + C>(^ 



(BIO) 



If we now use the expansion ( |3.7| ) of fi c , we obtain the result ( |3.9| ) presented in the main 
text. 



2. Cancellation of logarithms 

To understand the origins of the logarithms in the preceding analysis, consider a straight, 
naive, perturbative expansion in p. Treating the —pip*ip term of the Lagrangian as a per- 
turbation, the logarithms then arise from the diagrams of Fig. [TJ|. Each diagram should be 
understood as evaluated at fixed x, with effective chemical potential /i = p — V(x), and then 
the result of the diagram integrated over x. p is treated perturbatively, while V(x) is not. 
The imaginary-time propagators in this perturbation theory, derived from the action ( [2. 3D , 
are 

Go(po,p) = -. p } ■ - , , (Bll) 

where we now set U — 1 for convenience. The logarithms are produced by the infrared 
behavior, near the center of the trap, of the loops drawn large in the figure. Specifically, it 
is the po=0, p— >0, x— >0 behavior of these diagrams, where p is the loop momentum of those 
loops. The small loops are UV dominated and so, to this order in the expansion in a, are 
insensitive to x near the center of the trap. The infrared divergences due to the large loops 
then produce the same common factor for all diagrams: 



J d 3 x J d 3 p [G(0,p)] 3 oc J = log divergent, (B12) 



where we've introduced the 6-dimensional phase-space vector 




(B13) 



The cancellation of these logarithms at the phase transition occurs because, at the phase 
transition, the inverse susceptibility vanishes at the center of the trap. This condition 
is shown diagrammatically in Fig. 0, which implies that the logarithms generated by the 
diagrams of Fig. Ill cancel each other at the order of a under consideration. 
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FIG. 14. Diagrams producing the infrared logarithm in the small /2 expansion of N. Each diagram 
should be understood as being evaluated with an effective chemical potential fi = p, — V(x) with /2 treated 
perturbatively. The dots represent the 2-point vertex coming from treating the —p,ij)*yi term in the La- 
grangian perturbatively. The crosses are as in Fig. |[ 

APPENDIX C: NUMERICAL RESULTS FOR SUMS 

The following sums were computed numerically using iterative application of the Euler- 
MacLaurin formula. 

£ ,3/2,3/2,1- i ;n/2 ^ 2.416 942 200 (CI) 



z 3/2 j3 /2( i+J )l/2 
i3/2j3/2 



(j + ,U/2 _ ,1/2 _ ,1/2 
E 1 Jj , 3 /2,/2 - - 8 - 215 157 561 (C2) 



^(ijfc^i + fcj^ + fcj^+j + fcji/a " 2 - 2H1 (C3) 

E ( i<7 -)3/2 fc l/2 ((« + j + jfc)l/2 ~~ fcV^ " " 16 - 70 (C4) 



APPENDIX D: DERIVATION OF THE SUNSET DIAGRAM 



Using the methods outlined in Appendix A, we will now reduce the sunset diagram ( |4.21 ) 
to the integral representation ( |4.22| ) in terms of distribution functions. Starting with 

2\mJ Tqkl {tl +Ui){tq Q +u q ){tko + ouk) 
and rewriting the Kronecker 5 function in ( |D1|) as in Appendix A, we obtain 

1 /87ra\ 2 r@ _l e~ ial ° _r e +iaqo t e +iak ° 

j p Tq iq + u q : 

Using the frequency sums (|A4|) and then performing the a integration, 

1 /8vra\ 2 f , . . . , e P"i - e P("«+"k) 
/ niuJi) niujq) n(u k ) 

/qkl LOl -Ug- LO k 



nf )(0) = -i(^) 2 /^^f ^(2^)(.- q - k ). (D2) 

v 2 V m J Jo Tl tl + ujpJQ tqo + ujqJK tk + uJ k 



< ] (0) = -- — / n(a*) n(u; fl ) n(«*) (2vr)^(l - q - k) . (D3) 



The integrand is well behaved at u>i = u q + uj k (except for the infrared divergence where I, 
q, and k all go to zero, which is dealt with in the main text). However, as in Appendix A, it 
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is convenient to introduce a spurious principal part prescription at this stage. Then, using 
and permuting integration variables, one arrives at ( 4.22j ). 
The last term in (|4.22 ), involving just one n, is proportional to 



P.P. 

qkl LOi - LO q 



^L(2vr)^)(l- q 



k). 



(D4) 



The qk part of this integration is just a special case of ( A 1 1|) with the momentum labels 
changed and p set to zero. As described in Appendix A, it therefore gives zero contribution 
in dimensional regularization for d = 3. 
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